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Abstract 

Brooding over bosons, wave packets and Bose - Einstein correlations, we present a generic 
quantum mechanical system that contains arbitrary number of bosons characterized by wave- 
packets and that can undergo a Bose-Einstein condensation either by cooling, or increasing the 
number density of bosons, or by increasing the overlap of the multi-boson wave-packet states, 
achieved by changing the size of the single-particle wave-packets. We show that the n-particle 
correlations may mimic coherent or chaotic behavior for certain limiting wave-packet sizes. Ef- 
fects of complete 7i-particle symmetrization are included. The resulting weights which fluctuate 
between 1 and n! are summed up with the help of a formal analogy between the considered 
wave-packet system and an already explored multi-boson plane-wave system. We solve the 
model analytically in the highly condensed and in the rare gas limiting cases, numerically in the 
intermediate cases. The relevance of the model to multi-pion production in high energy heavy 
ion physics as well as to the Bose-Einstein condensation of atomic vapors is discussed. As a 
by-product, a new class of probability distribution functions is obtained and the critical density 
for the onset of pion-lasing is derived. The multiplicity dependence of single-particle momentum 
distributions are predicted that could be utilized in future event-by-event measurements. 



1 Introduction 

The study of the statistical properties of quantum systems has a long history with important recent 
developments. In high energy physics, quantum statistical correlations are studied in order to 
infer the space-time dimensions of the elementary particle reactions. The measured characteristic 
length scales are on the 10~^^ m scale, the time scales are of the order of 10^^^ sec. In high 
energy heavy ion collisions hundreds of bosons are created in the present CERN SPS reactions when 
Pb + Pb reactions are measured at 160 AGeV laboratory bombarding energy. At the planned RHIC 
accelerator, thousands of pions could be produced in a unit rapidity interval. If the number of bosons 
in a unit value of phase-space is large enough the bosons may condense into the same quantum state 
and a pion laser could be created |]l]. Similarly to this process, when a large number of bosonic 
atoms are collected in a magnetic trap and cooled down to increase their density in phase-space, the 
bosonic nature of the atoms reveals itself in the formation of a Bose-Einstein condensate which 
can be considered as a macroscopic quantum state. This condensation mechanism may provide the 
key to the formation of atomic lasers in condensed matter physics. 

Thus it is essential to understand how the quantum statistical properties of a dilute bosonic gas 
change if the density of the quanta increases. In terms of multi-particle wave packet states this 
corresponds to the increase of the overlap of the single-particle wave-packets and this is indeed the 
key mechanism in our picture to provide the Bose-Einstein condensation effects. Due to the fact 
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that the quanta in high energy physics are created by well localized sources, we would like to avoid 
the plane-wave approximation and utilize wave-packets instead. In case of the Bose condensation of 
atomic vapors, the presence of the magnetic trap prohibits the application of the outgoing plane- 
waves. We hope that the wave packet formalism, presented herewith, may have certain applications 
in that area of condensed matter physics as well, due to the finiteness of the wave-packet sources. 
However, we will not explore this possibility in details herewith, since our main motivation for 
considering the problem of quantum statistical correlations for many hundreds of overlapping wave- 
packets stems from the experimental program in high energy heavy ion physics, which itself has 
been motivated by a search for a new phase of the strongly interacting matter. This new phase is 
referred to as the Quark-Gluon Plasma (QGP), where the quark and gluon constituents of strongly 
interacting particles are expected to be liberated from the confining forces of the strong interaction. 
See refs. Q for recent reviews on the experimental and theoretical developments in the quest for 
the QGP. 

Bose - Einstein correlations have already been investigated in a large number of papers for various 
reactions of high energy physics (see for example the list in refs. 0). In particular, some models 
assume that pions are created at a given point in space with prescribed momenta that can only 
be approximately valid in a quantal description. Such models correspond to certain semi-classical 
Monte Carlo simulations of high energy reactions, e.g. [|[ |lC, 11, 1^. The Bose-Einstein correlations 



are switched on only after the particle production is completed in refs. [p| , |lO| , |13| with the help 
of a weighting method that utilizes the Wigner-function formalism In the latter, the 

production of the pions is given in space-time as well as in momentum space, since the hadronic string 
model associates a momentum space fragmentation with a space-time fragmentation, see ref. [ p^ . 
There are different approaches incorporating both quantum mechanics and source models for pions, 
such as the Wigner function formalism, developed by Pratt |14[ , and its different generalizations, 
for example the relativistic Wigner function formalism of ref. [ |15[ . 

In the present paper we will utilize wave packets in order to resolve the clear contradiction of 
the classical assumptions in some models of Bose-Einstein correlations with the tenets of quantum 
mechanics. Although certain properties of the wave-packets have been already utilized in previ- 
ous publications |T^, these were usually applied on a phenomenological level, utilizing e.g. 
minimal wave packets on the single particle level to mimic the effects of the uncertainty relations. 
However, in our case, we shall consider the problems arising from the overlap between different 
wave-packets, we shall discuss, how a Poisson multiplicity distribution is gradually transformed to 
a multiplicity distribution of a Bose-Einstein condensate for a suitable choice of the wave packet 
size. {Note that this type of Bose-Einstein condensation is different from the text-book example of 
Bose-Einstein condensation in statistical mechanics, where the zero-momentum plane-wave modes 
are populated with a macroscopic amount of quanta. Here the wave-packet state with zero mean 
momentum will be populated with a very large number of particles. Due to the optical coherence of 
the bosons in the zero mean momentum wave-packet state, this system is sometimes referred to as 
a laser, c.f. the term "pion-laser" introduced by S. Pratt in ref. 

We shall also see how the singe-particle spectrum and the two-particle correlation functions shall be 
drastically changed if the overlap of the multi-particle wave-packets reaches the critical size necessary 
to create the Bose-Einstein condensate. 

In high energy heavy ion reactions the system at break-up is considered as a pion source, which in 
general can only be described by a density matrix. One may describe a density operator by a weighted 
sum of projection operators. We will follow this method, and specify the various pion sources in 
two steps. In section ^ we summarize the generic formulae of the density matrix formalism and 
summarize the definitions of the physical quantities under investigation. At this level, projectors to n- 
particle wave-packet states shall be introduced. These shall be constructed from single-particle wave- 
packet states in subsection 2.1. In subsection 2.2 we specify a density matrix that includes induced 
emission on the level of multi-particle wave-packets. We reformulate the multi-particle wave-packet 
density matrix model in terms of new variables in such a manner that an effective multi-particle 
plane-wave model is obtained, with source parameters explicitly depending on the wave-packet size. 
We utilize of the ring-algebra of multi-particle symmetrization - discovered by S. Pratt [Q and 
described in details in ref. 0] - to reduce the effective plane-wave model to a set of recurrence 
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relations in subsection 4.1. We solve these recurrences analytically: A subset of the recurrence 
equations is solved analytically in the case of arbitrary boson density in section ^ . Complete analytic 
solutions are obtained for the generating function of the multiplicity distribution. Analytic results 
for the single particle inclusive momentum distributions and the two-particle inclusive correlation 



functions are presented in subsections 5.1 and 5.2 for a rare gas limiting case and a highly condensed 



Bose-gas limit, respectively. The onset of the wave-packet version of Bose-Einstein condensation is 
studied analytically in section |5.3| . The analytical results presented in these sections are new results 
not only for the multi-particle wave-packet systems but also for the multi-particle symmetrization 
of plane-wave systems. In section ^ we present numerically evaluated single-particle spectra and 
two-particle correlation functions for the generic case. Finally we summarize and conclude. 

2 Formulae for a generic case 

The density matrix of a generic quantum mechanical system is prescribed as a sum of density 
matrixes with different number of bosons, 

oo 

P^^PnPn, (1) 
n=0 

the density matrix of the whole system is normalized as 

Tr (p) = 1, (2) 

similarly to the density matrix characterizing systems with a fixed particle number n, which satisfy 
the normalization condition 

Tr (/3„) = 1, (3) 

from which it foUows that 

oo 
n=0 

The multiplicity distribution is prescribed in general by the set of {pn}^o- Later on we shall utilize 
the Poissonian multiplicity distribution for the case when multi-particle symmetrization effects can 
be neglected. 

The density matrix describing a system with a fixed number of bosons is given by 



/5n = y (iQi...da„ p„(ai, a„) |q!i, Q!„)(ai, a„|, (5) 

where the states |ai, a„) denote n-particle wave-packet boson states, which shall be prescribed in 
more details in the next subsection. The parameters ai characterize a single-particle wave-packet, 
and the multi-particle wave-packet states are properly normalized to unity: 

(ai, Q;„|ai, q;„) = 1. (6) 

Thus the normalization condition for the n-particle density matrix can alternatively be written as 

J dai...dan Pn{ai, ...,an) ^ 1- (7) 

The expectation value of an observable, represented by an operator O is given by 

<6>=Tr(/56) (8) 

The i-particle inclusive number distributions are given by 

iV,(ki,...,kO = Tr (pat(k^)...„t(k^)a(k^)...a(k,)) (9) 
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which is normahzed to the i-th factorial moment of the momentum distribution as 

d^ki ...(fk,N,{ki, ...,k,) = {n{n~ l)...{n - i + 1)). (10) 

These inclusive number distributions can be built up from the exclusive number distributions of 
A^™(ki, k„i) for all m = « + 1, ... describing the i-particle number distribution for fixed m 
multiplicity events as 

oo 

iV,(ki, k,) = K^i^^lki, kO (11) 

m—i 

where 

iV^"^(ki,...,kO =Tr (p™at(ki)...at(k,)a(k,)...a(ki)) (12) 

and 

J d^ki...d^k,N^"'\ku...,k,) = m...{m-i + l) (13) 
The i-particle probability distributions are defined as 

P^(ki,...,k,) = -^4 — TT7A^»(ki,...,kO (14) 

{n[n — l)...[n — i + 1)) 

which are normalized to unity. For a fixed multiplicity of m the i-particle probability distributions 
are defined as 

P^)(ki,...,k,) = — -4 ^7V('")(ki,...,k,) (15) 

m{m — l)...(m — I + 1) 

which are also normalized to unity. The i-particle inclusive correlation functions are defined as 

a(ki,...,kO= Pf'-'^^^ (16) 

riA(k,) 

while for a fixed multiplicity m the i-particle exclusive correlation functions are defined as 

ci"'^(ki,...,kO^ 1^1. -.J^O (^7) 

ri^^^(k.) 



Let us note, that another definition of the i-particle (exclusive/inclusive) correlation functions has 
been discussed recently in ref. Q 
respectively. These are defined as 



been discussed recently in ref. M, which will be denoted by Cf'^"^\ki, ...,ki) and Cj^iki, ...,ki), 



Cf-(™'(k„...,kO = (18) 

f[Nt\k,) 

Cf(ki,...,kO = ''f^-^^l (19) 

ri^i(k,) 

and they are related to the probability correlations by a normalization factor, 

cf^")(k„...,kO . (20) 

m 

cf(k„...,k.) . ^^^^-']-^:^-' + '^^ aik,^...,k.). (21) 

(m) 
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These relations indicate that for large values of m the difference between the exclusive correlation 
functions of eqs. (18) and ( p^ is small, but the definition may matter much more when the inclusive 



correlation functions are evaluated from eqs. (^9|) or (|T^. 
The one- and two pion distributions are calculated as 

7V2(ki,k2) - Tr (/5at(ki)at(k2)a(k2)a(ki)) 

iV2(kl,k2) = ^p„iV2("^(ki,k2) 



7vf)(ki,k2) = TV (/5„at(ki)at(k2)a(k2)a(ki)) (22) 



Similarly, 



iVi(ki) = ^p„7V(")(ki) 

11 

7V("'(ki) - Tr (p„at(ki)a(ki)) (23) 
We also have the following general relations 

Jd'k2Nt\k,M) - (n-l)7v(")(ki) 

y"d3kiiv}")(ki) = n 
d^kid^k2N2{ki,k2) = <n{n-l)> 

J d^kiNiiki) = <n> (24) 
The two-particle inclusive Bose-Einstein correlation function is then given by the expression 

r(k k)-r(k k)- ^2(ki,k2) _ P2(ki,k2) 

< n{n- 1) > iVi(ki)A^i(k2) Pi(ki)F(k2) 

The general considerations presented above provide the framework for including the effects arising 
from multi-particle Bose-Einstein symmetrization using wave packets. Before specifying how the n- 
particle density matrixes /5„ are related to the n-particle wave-packets, let us first introduce single- 
particle wave-packet states. The simplest case, when a single particle is emitted from a single 



source, shall be specified in subsection 2.1 and shall be extended for many particle wave-packet 



states emitted from many sources in subsection 2.2. We end this subsection by specifying a density 



matrix for which one can overcome analytically the difficulty related to the over-completeness of 
single-particle wave-packet states which shall complicate the normalization of multi-particle wave- 
packet states. The higher-order symmetrizations are reduced to an equivalent plane-wave problem 
in section ^, that is solved analytically in section |^ in certain limiting cases. Numerical results are 
shown in the subsequent section || to illustrate the new effects in the generic case. 

2.1 Single-particle wave-packets 

The building block for the definition of states will be a single particle wave packet, which is described 
as follows. The operator a^(x) creates a boson (e.g. a pion) at point x. This creation operator at(x) 
can be decomposed as 

-—3- exp(ipx)at(p). (26) 

The creation and annihilation operators in space and in momentum space obey the standard com- 
mutation relations 

[a(x), at(x')]=^(x-x') (27) 
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[a(p), at(p')]=^(p_p') (28) 

Let us define a wave packet creation operator as a linear combination of the creation operators in 
space, 

= / exp[-(x-e.)'^V2] exphi^.x] a^^). (29) 

Here ai — (f;, tt;, at) refers to the parameters of the wave packet: the center in space, in momentum 
space and the width in momentum space, respectively. (The index i is introduced to distinguish 
single-particle wave-packets with different parameters). The wave-packet created by this operator 
is given by 

I«.>=«1|0), (30) 

which is normalized properly as 

{a.,\a,) = 1. (31) 

Although these wave packages are normalized to unity they are not orthogonal; in fact they form 
an over-complete set pl[ . The commutator for the wave-packet creation and annihilation operators 
is given by the overlap as: 

as can be verified directly from the definition of the wave-packet creation and annihilation operators. 
The overlap between two wave-packets is given explicitly by 

(3/2) 



[a^,a]] ^ {a,\aj), (32) 



' 2(a,^ + a|) ) 

exp I i{TT^ - TTj) I ' '"2 ', 2"' I I (33) 





The coordinate space representation of a wave packet is given by 



exp[-^ ^ iTTx] (34) 



Expansion ( p6|) can be used to obtain the momentum space representation as 

(Pl«> = -4v^ cxph^P-^ - t^P - n)] (35) 

Up to this point all quantities were defined at a given time tg, which may be considered as break-up 
time. The notation la) = \a, to) refiects this fact. Assuming that after the break-up the wave packets 
move without interaction, the time evolution of the wave packet state is given as 

\a;t) = e'^(*-*°) \a) = J d'p{p\a) exp[zc.(p)(< - <o)] a^(p) |0) (36) 
We introduce the notation 

w{a, p) = {p\a; t) = ^(p, tt) exp[iw(p)(i - io) - i(7r - p)C] (37) 

with 

1 r (P — Tt)^ n , , 

g(p,^) = — ^exp[- ^^ J ]. (38) 

The pure state in eq. ( ^^ describes a wave packet with mean momentum tt^, mean position x^ at 
time Ti, with a spread of cr = Ap in momentum space and a corresponding Ux = Ai? = Ti/a spread 
in configuration space. Although the location and the momentum of a pion cannot be specified 
simultaneously with an arbitrary precision, the values of the parameters 7r,^,cr are not subject to 
any quantum mechanical restriction. For example, these values of the wave packet parameters can 
be identified with the classical pion production points in the phase space, as calculated from a Monte 
Carlo or a hydrodynamical model. 
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2.2 Building multi-particle wave-packet states 

The basic building block for the following considerations is the operator al{t) which is defined by 

alit) = a| eMMp)it - h)] (39) 

where the index i refers to the difi^erent possible values of the center of the wave packets in coordinate 
and momentum space. For simplicity we may assume that all the wave packets are emitted at the 
same instant. Thus we indicate (TTi, ^i, crj, to) by a^. Later on when we evaluate the recurrence 
relations we shall also assume that each wave packet has the same size, a = ai. 
Let us assume that there are M boson sources in the system. If one boson can be emitted from any 
of these sources then the state vector can be written as 

M 

|{a}; 1; i) = |ai, c^m; 1; t) a ^ al(<) |0). (40) 

i=l 

This state is a linear combination of one boson states and is not normalized. The normalization 
constant is to be determined from the condition ( {aj}; 1; t\ {a.^}; 1; t) = 1. 
The n-boson states can in turn be given as 

M 

\{a, }; n; t) = \ai,...,aM; n; t) cx a|(t)]" |0). (41) 

i=i 

In the usual treatment M ~ n. For example, if M = n = 2, the two-particle state is specified as 

|ai,a2; 2 t) oc [ a\{t) + al{t)] [ a\{t) + al{t)] |0). (42) 



When using this type of states one usually introduces a phase-averaging method to get back the 
Bose -Einstein enhancement of pions at small relative momenta . A different possible n-boson 
state is given by 

L 

{MJ; n; t) = a^; Mi,...,Afi; n; i) cx [][ "I (*)]*'' |0). (43) 

where 

Mi+M2 + ... + ML = n, (44) 

and Ui Uj for i j- 

We restrict our investigation to this latter case only. In eq. (^) there are exactly n creation operators 
acting one after the other on the vacuum, and the identical creation operators are grouped together. 
The most general case is that similar wave-packets can be created repeatedly and that it is not 
necessary to create two similar wave-packets one after the other, but a different wave packet can be 
created in between. Because of this reason, the general form of eq. (|3|) can be written as 

n 

|ai,...,a„;t) (X n"l(i)|0) (45) 

i=l 

where similar values are allowed for arbitrary subsets of {a^}"^]^. As a special case, it is possible 

that all bosons are created in the same wave-packet state. 

The norm of the state defined in eq. ( p5| ) can be calculated using the fact that 

n n 

(0|a(pi) ... a(p„) at(p:j ... at(p;)|0) ^ ^ [] WpO, "^(p^J] - ^11 '^(p^ - p'.J (46) 

cr(") i=l ct(") i=l 
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where the summation over cr^"'' indicates the summation over the permutations of the first n positive 
entire numbers and o-. denotes the permuted value of i in a given permutation. ( Note that the 
subscript ^. which stands for the permuted value of i should not be confused with the normal-size 
(Ti, which stands for the width of the i-th wave packet in the momentum space.) Similarly to the 
plane wave case, one has 

n n 

(0| ai ... a„ ai ... a\ |0) = ^J] = EH ("^1"--) (^7) 

o-(") i=l <t(") i=l 

The n boson states, normalized to unity, are thus given as 

1 



Qfi, ... , an) = I ai, a„; t) 



aiit) ... a\{t)\0). (48) 



Let us introduce the notation 

7, J = {a,\aa^) (49) 

which can be further simplified if the width of the wave-packets are identical, ai = Uj. This 
simplifying condition is in principle not necessary, but we shall apply it in the forthcoming because 
this assumption will simplify certain recurrence relations. If all the wave packets have the same 
width, one has 



li,3 



J (fpw*{p,ai)w{p,aj) 

exp( ^-^)exp( ) X 

exp(i7Ti(,j - i-Kjii) exp(i(7r,; + 7rj)(^i - Ci)/2) (50) 



Note, that 7i_i = 1 and 7ij- = 7j'i- The summation over the permutations cr*-"^ can be decomposed 
as summing over transpositions. Since any permutation in ct^"^ can be built up from the product of 
at most n — 1 transpositions (i, j), the summation over the permutations can be written as a sum 
over the partial sums where a partial sum contains all the terms which belong to a class of a given k 
with < fc < n — 1, where k is the minimal number of transpositions necessary to build up a given 
permutation: 

n 

(ij) (»,i);0-,fc) (»j);(fc,0 

For n = 2 and n — i the normalization factors are given, respectively, by the inverse square root of 
the quantities 

1+ I 71,2 ? (52) 

1 + I 71,2 P + I 72,3 P + I 71,3 P + 71,2 72,3 73,1 + 7l,3 73,2 72,1 (53) 

In the next section we shall introduce a new type of density matrixes using the above defined set of 
normalized wave-packet states. Later on we shall provide the analytical solutions to these models. 



3 A New Type of Solvable Wave-Packet Density Matrix 



In the formulae displayed up to now the normalization factor always appears as e.g. in eq. (51). For 
a reaction like Pb + Pb at SPS in a single central collision event one has approximately 500 - 600 
TT^ . Thus for this case we have the the sum over all possible permutation of the 600 pions, and 
each term in the sum is the product of 600 functions, depending on 1200 vector parameters. And 
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this expression is in the denominator. Later on, this normahzation factor will be multiplied with the 
distribution function of these parameters, and than one should integrate over these parameters. It is 
clear, that it is practically impossible to perform this task. A very interesting numerical simulation 
has been reported recently about a very reasonable numerical approximation to evaluate the n-boson 
symmetrization effects with the help of a Monte-Carlo algorithm However, the efficiency of such 
algorithms decreases as l/{nl) thus it is practically impossible at the moment to go beyond 5-th or 
6-th order explicit symmetrizations. 

There is one special density matrix, however, for which one can overcome such a difficulty even in 
an explicit analytical manner. Namely, if one assumes, that we have a system, in which the emission 
probability of a boson is increased if there is an other emission in the vicinity. This would be an 
effect similar to the induced emission. Such an n-boson density matrix, eqs. may have the 

form: 

p„(ai,...,a„) = -^^^ [ J|pi(aO ) ( IT ) (^'^^ 



/ fe=l 



The coefficient of proportionality, M{n), can be determined from the condition that the density 
matrix is normalized to unity, 

/ n n „ n 

Wda^pi{ai)Y^ n " XI / Wd^^PM^){'^k\a<yk) (55) 

i=l cr(") k=l (t(") i=l 

The density matrix given by eq. (|4|) corresponds to the expectation that the creation of a boson has 
a larger probability in a state, which is already filled by another boson. The above model not only 
describes such a source but it makes possible to continue the calculation as well, given that the last 
term of the density matrix cancels the normalization factor of the overlapping wave-packet states. 
The induced emission, that is implicitly built in to the above definition of this density matrix, 
can be made much more transparent, if one evaluates the ratio pn{cti, Q;„)/[n"=i Pi(Q^i)] for two 
special cases: one when each of the n particles are emitted with the same wave-packet parameter ai 
(maximal overlap) and the other, when the overlap between the wave packets of any pair from the 
n can be negligible. If n particles are emitted in the same wave-packet state, one has 

p„(ai, ...,ai) n\ 

(56) 



[pi(ai)]" N{ny 
while if the overlap between any of the wave packets is negligible, we obtain 

p„(ai, a„) 1 



N{n) 



(no overlap). (57) 



In general, the overlap of the wave-packets determines the magnitude of the enhancement of the 
density matrix: 



p„(ai, ...,q;„) fe=i 



AA(n) 



(overlap). (58) 



Thus the density matrix given in eq. (|5j) describes a quantum-mechanical wave-packet system with 
induced emission, and the amount of the induced emission is controlled by the overlap of the n 
wave-packets, yielding a weight in the range of [l,n!]. Although it is very difficult numerically to 
operate with such a wildly fluctuating weight, we shall show that the special form of our density 
matrix yields a set of recurrence relations that can be evaluated numerically in an efficient manner. 
The essence of the method is the reduction of the problem to the already discovered "ring" - algebra 
of permanents for plane- wave outgoing states j^. 
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We have a larger freedom to choose the form of pi (a) , the density matrix describing the distributions 
of the parameters of wave packets. However, one can perform the calculations analytically when 
choosing Gaussian forms. For simplicity, we do not discuss the fluctuations of the wave-packet 
sizes in the forthcoming, although a fluctuating wave-packet size with Gaussian random distribution 
would lead to a similar mathematical structure of the solution. 
Thus for the single-particle density matrix we assume the following form: 

Pi{a) = p^^iO Pp{Tr) S{t - to) 

= (2.Jt)3/2 -M-^Vi2mT)) (59) 

Note, that this choice corresponds to a non-relativistic, non-expanding static source at rest in the 
frame where the calculations are performed. 

This completes the specification of the model. // one is interested only in the numerical results, 
we recommend to jump directly to Section where multiplicity distributions, particle spectra and 
correlations are plotted from a numerical evaluation of the model. For those more theoretically 
interested, we present in the next sections the analytic solution of this model. We may comment 
that these solutions are not easily obtainable, however, the algebraic beauty of the structure of the 
equations well compensates for the difficulty. 

4 Algebraic Evaluation of the Model 

In this section we make algebraic manipulations that are necessary to reduce the multi-particle 

wave-packet problem to an analytically solvable plane-wave problem. 

The n particle distribution for the system containing n bosons is given by the expression: 



^n"-* (ki,k2, ...,k„) = —:-—/ ]^ da™yo(a„)^J| w*(ki, a,;) w(ki,a^J 

^'"''■^ m=l i=l 



(60) 



Since the indices run through all possible values, we can exchange the index of k and a in the second 
w. 

iV^") (ki,k2, ...,k„) = — :— / ]^ damp(am)^ J|-u;*(ki,Q;i)«7(k^^,Q;i) (61) 

m=l ct(") i=l 

One more rewriting becomes possible by using the identity (11"=! ^0(11^=1 ^j) — Yli=ii^i^i)- 

1 " /■ 

iV^")(ki,k2,...,k„) = —^yT] da,w*{k,,a,)p{a,)w{K,,a,) (62) 

Let us introduce the auxiliary quantity 

p(kj,kj)= daip{ai)w*{'ki,ai)w{kj,ai) (63) 



where the overline and the two arguments distinguish this auxiliary quantity from the single-particle 

density matrix p(a), introduced earlier. 

Using these notations Eq. (p3) can be rewritten as: 



7V(")(ki,k2,...,k„) = ^-^^J]p(k„k,) (64) 

^ > i=l 
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With this notation the "source function" will be defined as 

Sn{ki,...,k„)^Y.I[p(^^'^-'^ (65) 

and further: 

iV(")(ki,k2,...,k„) = -^5„(ki,...,k„) (66) 

Using Eq. (^) and the definitions (|l^,|l^), the one and the two particle distributions are cast in 
the form 

7Vl"^(ki) = J []dk;5„(ki,...,k„) (67) 

4"\k,,k2) = "^r^ I ndkj5„(ki,...,k„) (68) 
From this it is clear that the constant of normalization M{n) can be expressed as 

/ n 
[]dkz5„(ki,...,k„). (69) 
1=1 

At this point one can realize, that the above equations have algebraic structure similar to those 
discussed in Refs. [Q, ^ , thus we can use their method to perform the summation over the permu- 
tations. 

Up to this point the expressions in this subsection refer to events with fixed pion multiplicity. 
For an average over many events, one has to treat the multiplicity distributions too. 
Following Refs. ^ we introduce the auxiliary parameter, riQ, which correspond to the mean 
multiplicity of a source containing classical particles (not symmetrized system). If the system is 
characterized by a Boltzmann distribution, then this no is a function of chemical potential, the 
temperature and the volume (characterized by the radius, R. ) 

Let us assume a Poissonian multiplicity distribution for the case when the Bose-Einstein effects are 
switched off (denoted hy pn^), i.e. let us assume that 

pW = ^exp(-no), (70) 
n! 

Let us define 

Gi(p,g) = ^p(p,q) (71) 

G2{p.q) = p(p,pi)dpip(pi,q) (72) 

Gi{p,q) = / P(P,Pi)c^PiP(Pi,P2)---p(Pi-2,Pj-i)dPi-ip(P*-i:q) (73) 

With these definitions we find that the case of multi-particle wave packets can be considered as 
a formally equivalent plane -wave system due to the complete formal analogy between the above 
equations and those of ref. . 

The above mentioned problem was shown to be solvable with the help of the so-called ring algebra 
discovered first by S. Pratt ||l|. Utilizing that algebra, a set of recurrence relations were obtained 
that contained some effective parameters of the T matrix which were assumed to have a Gaussian 
form. In our case, the width of the single-pion wave packet and the width of the distributions of the 
wave-packet centers are taken into account explicitly, thus a wave-packet size appears as one of the 
parameters of the recurrence given below. We shall also present certain analytic solutions to these 
recurrence relations and investigate them numerically as well. 
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4.1 Recurrence relations 

Having performed the reduction of the multi-particle wave-packet problem to the multi-particle 
plane-wave problem, we now can utilize the reduction of the multi-particle plane-wave problem to a 
set of recurrence relations that were given first in rcfs. 

The multiplicity distribution (with the inclusion of the symmetrization effects), the one- and two- 
particle distribution for a system containing n pions can be expressed with the help of the three 
definitions that relate these observables to the elements of the recurrence relations, four recurrence 
relations and four initial conditions (one for each recurrence relation). Two further equations are 
necessary to define certain auxiliary variables that are related to the recurrence relations and the 
observables. Thus one ends up a set of 14 equations, as given below. The observables are defined as 

Pn = (74) 



fe=0 



1 " 

p(»)(k) = Vc^„_,G,(k,k), (75) 



2—1 



- n l—l 

Pt\ki,k2) = T^Y^Y^ ^n-l [G™(ki,ki)Gz-™(k2,k2) +G™(ki,k2)Gi-™(k2,ki)] , 

n(n-l)a;„^;^^ 

(76) 

in order to relate these definitions to the recurrence relations, one introduces two auxiliary quantities 
as 

G„(p,q) = < /i„exp(-a„(p2 +q2) + g„pq), (77) 

G„ = -fd^pGUp.p), =K^(—^ y . (78) 

n J n \ 2an- On/ 

The physical interpretation of the quantities LUn and G„ shall be discussed later. The recurrence 
relations correspond to the solution of the ring-algebra |l| are given as 



1 " 

-YlClUJn-l (79) 
7^3/2 

^ (80) 

(ai -|-a„) ' 

fln+i = ai - ^\ r (81) 

4(ai -f a„) 

9n+l = -TTf 7 (82) 

2[ai + an) 



The initial conditions or the starting elements of the recurrence are 



^0 = 1, (83) 
1 

91 = Riff, (86) 



which yield the following value for Gi : 



Gi = no. (87) 
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In the initial conditions eq. ( p4[j86| ) the foUowing notation is introduced: 

(t|i = (T^ + 2mT, 

TflT 

^eff ^ ^ ^ T^- (89) 

or using the spatial spread, ax, and the effective temperature T^ff, that characterizes the system 
before symmetrization effects are taken into account, one can write 



Teff - T + a'/{2m)^T+l/{2mai), (90) 

Riff = ^' + 4^- ^^^^ 

2 Teff 



Observe that the above recurrence relations correspond to the pion laser model of S. Pratt when a 
replacement R — s- Reff and T Teff is performed . The parameters of the T matrix in S. Pratt's 
model are thus interpreted here as effective parameters and expressed with the help of the radius of 
the source of the wave-packets, the temperature of the source of the wave-packets and the width of 
the wave-packets in our multi-particle wave-packet model. 

Thus the reduction of the multi-particle wave-packet model to an equivalent plane-wave model is 
completed. These recurrence relations were studied numerically but analytic solutions to these 
equations were not given before. Although numerical investigation of the presented recurrence 
relations were reported e.g. in refs. ^, to apply these kind of numerical investigations for 
multiplicities of few 1000 high numerical precision and special technique is needed to treat the 
factorials of large numbers appearing in the algorithm. We present some new numerical results in 
this multiplicity range in Section ^ 

Before proceeding to the solution of the above recurrences, let us note, that in eqs. ( p^[76| ) the 
multiplicity distribution and the exclusive single-particle and two-particle momentum distributions 
were given only, i.e. the particle spectra that are measured, which belong to special events with fixed 
multiplicity n. Such events occur with a probability of p„, i.e. they constitute a sub-set of all the 
events. To get the more readily observable inclusive momentum distributions ~ hence the inclusive 
correlation functions - one has to average over all the multiplicities with weight factors p„. To our 
best knowledge, such an averaging of the exclusive momentum distributions was not performed before 
in similar kind of models, neither analytically nor even numerically, due to numerical difficulties 
involved in evaluating the exclusive distributions for a large number of n. However, the following 
analytical considerations can be utilized to simplify this task. Let us consider the identities 

oo / n \ 
= I] II a«-m?>m (92) 




n=2 \m=2 r=l ) 

When inserting the expressions for the exclusive momentum distributions as given by eqs. ( |74| - 176| ) 
to the definition of the inclusive momentum distributions in terms of the exclusive ones, eqs. (|ll[-[l5[), 
exactly this type of products of infinite sums appear, in the form of the r.h.s. of the above relations. 
When rewriting these in terms of the l.h.s. of the identities, one of the factors shall be 'Y^=^Vn = 1 
that simplifies [|l] the result as 

oo 

7Vi(k) = ^G,(k,k), (94) 
i=i 



iV2(kl,k2) = K^G,(ki,ki) K^G,(k2,k2) +^G,(ki,k2) 



:G,(k2,ki) (95) 
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These simple and beautiful relations will be useful when evaluating the inclusive number distributions 
analytically as well as numerically. The probability distributions can be obtained from eqs. (p^). 
The importance of the proper normalization of the two-particle inclusive correlation functions was 
emphasized recently in ref. [Q. We are now in the position that a formula of general validity can 
be derived that involves no approximation and can be obtained without any reference to the details 
of the source model that determines the functions 6*^(141, k2). As long as the n-particle exclusive 
number-distributions have the form of eqs. - 16^), the probability distributions are given by the 
functional equations of eq. (^- ^) and the above eqs. - 19^) follow. In such cases it is possible 
to introduce the auxiliary quantity 

oo 

G(ki,k2) = ^G,(ki,k2). (96) 

(In the Wigner-function formalism, this quantity corresponds to the on-shell Fourier-transform of 
the emission function S'(a;,k).) With this notation, both kind of two-particle inclusive correlation 
functions, eqs. (|l^,|l^) can be evaluated as 

cfc,..) - ,98) 

(n(n-l))V G(ki,ki)G(k2,k2)/ 

Thus we find, in agreement with ref. [ po[ , that in general a non-trivial pre-factor appears before 
the two-particle inclusive momentum distribution G2(ki,k2), that can be transformed away if the 
correlation function is defined as the ratio of the number of counts, eq. (|l9|). Alternatively, this result 
can be interpreted in another manner: when measuring the inclusive correlation function as the ratio 
of the detection probabilities, eq. (|l^), an overall normalization constant has to be introduced that 
will be not arbitrary but will have to be equal to (n)'^ / {n(n—l)) , the inverse of the second normalized 
factorial moment of the multiplicity distribution. 

In the next sub-section we re-write these recurrence relations into dimensionless forms and solve 
analytically three of them completely and present a formal solution for the fourth recurrence too. 
We present the leading order complete solution of the model both in the rare gas and in the very 
dense Bose-gas limiting cases. 



5 Analytic Solutions 

When one assumes a Gaussian form for the wave-packets and for the distribution functions of 
the parameters of the wave-packets, the functional equations for G„-s are reduced to the set of 
recurrence relationships presented in the previous section. The recurrence equations for the exclusive 
distributions were known since their discovery in 1993 (except for the inclusion of wave-packet sizes 
into the definition of Reff and Tgff). Thus, a number of properties of these recurrences were 
explored but as far as we know, only numerically. In particular, it has been observed that at 
a critical density Uc the stimulated emission of bosons over-compensates for the decrease of the 
unsymmetrized pl"^ probabilities and a lasing effect or coherent behavior appears. The condensation 
is characterized by the divergence of p„ (symmetrized) probabilities with increasing values of n, by an 
appearance of a low-momentum peak in the single-particle spectrum and by a decreasing intercept 
of the correlation function: G(k, k) < 2. However, no analytical results were known about the 
spectrum and the correlations at the condensation point as well as in other limiting cases, only 
the exclusive correlation functions were evaluated numerically but the evaluation of the inclusive 
correlations was not performed. 

Note also, that until now uq was just the mean number of bosons before the symmetrization, and 
the parameters uq, Reff and T^ff were just interpreted in terms of the theoretical input values of 
the correlated system. It seemed that the actual single-particle spectra and correlation function can 
only be determined by numerically solving the recurrence relations. The invention of the recurrence 
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relation in ref. Q was already a huge step forward, since the number of steps that are required for the 
solution of the n-particle symmetrization in general increase as n\ cx exp[n(lnn— 1)] which increases 
faster than any polynomial - or a non-polynomially (NP) hard problem. The NP-hard problems are 
very difficult and inefficient to handle numerically for large values of n. S. Pratt reduced this difficult, 
NP-hard case to a set of recurrence relations, where the number of necessary steps to evaluate the 
observables increases only as slowly as and thus the problem is solvable even for large values of 
n ~ 1000 within a few minutes on the current computers. 

Here, we would like to present the first analytical solution for the multi-particle symmetrization 
of the wave-packets — no numerical evaluation of the recurrences will be necessary, and the time 
needed to compute the result will be independent of n. The solution of the problem is possible in 
the rare and the dense gas limiting cases, while in the general case the multiplicity distribution shall 
be given in terms of its combinants. 

Before presenting the analytical solution of the recurrence equations, let us re-formulate the recur- 
rences in eqs. (|o| - ^2|) for new, dimensionless quantities. Note that eq. (|79| ) is already referring to 
dimensionless quantities thus there is no need to reformulate it. 
Let us introduce the following dimensionless variables: 



A„ 
H„ 

Gn 



o-^(ai -I- a„) 

hn/hi = ;i„/(w|)3/2 



R 



eff'^T 



(99) 
(100) 
(101) 
(102) 



Note that the variables indicated by the upper-case An , Hn and G„ are essentially the dimensionless 
versions of the variables an, hn and gn- The quantity x corresponds to a dimensionless measure of 
the phase-space available for a single quanta. Extremely rare gas corresponds to the limit x oo 
while a very dense Bose-gas corresponds to the a; ^ limiting case. Note also that the dimensionless 
variable x should not be confused with the vector x indicating a position in space. 
The recurrence relations (pQ - p3) can be re-written for the dimensionless variables as 



Hn+l 
Gn+1 



A, 



x^J_ 

4 An' 



- Hnj-, 

— Gn 



(103) 
(104) 
(105) 



The initial conditions for these recurrences read as 



Ai = l + x, 
Hi - 1, 
Gi = X. 



(106) 
(107) 
(108) 



The dimensionless recurrences of eqs. (103-105) can be solved exactly with the help of the auxiliary 
quantity 



Yn — IT^^' 



i=l 



since the solution of eqs. (104-105) is given as 

Hn+l = 
Gn+l = 



y(3/2; 

x^ 



(109) 

(110) 
(111) 
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where the initial conditions of eqs. (107-108) are already taken into account. The remaining recur- 
rence relation, eq. (103) is easy to solve if one rewrites this as 



y„+i = (1 + x)Y^ - —Yn-u 



(112) 



n+1 



(l-fa;-7) 



n-l 



that can be further re-formulated for an arbitrary value of the parameter 7 7^ 1 + x as 

^ 

4(1 + a: -7)" 

and this equations becomes solvable if the value of the parameter 7 is chosen such that 

™2 



7 



This equation has the following two roots: 



4(1 + a: -7)' 



(113) 



(114) 



7± 



(1 + a; ± VI + 22;) 



(115) 



Note that these roots satisfy the following useful algebraic relations: 



7+ +7- 
7+7- 

^(1/2) _ ^(1/2) 



1 + 

1. 



(116) 
(117) 
(118) 



One may choose any of 7+ or 7_ to solve the third recurrence equation by re-scaling eq. (112) and 
introducing 



Vt = Yr,+ i - -I^Yn 



(119) 



One finds that 



2/1 



7±y^-i, 



(120) 
(121) 



thus the solution for the third recurrence relation can be written as 



(122) 



which yields the solution for Yn as 



Yn 



7± - 7t 



(123) 



that can be substituted back to eqs. (110-111) to get an explicit solution. Then the variables that 
have dimensions can be also calculated easily, and their solution in terms of x = R^j^a^ is especially 
simple if the variables defined in cq. (115) are used: 

7+ - 7- 7+ + 7- 



hn 



9n 



20-2 

1 



72.13/2 



7+ - 7- 



7?-7- 



(3/2) 



CJj' 



(124) 
(125) 
(126) 
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and these solutions can be utilized to evaluate the quantity C„ from the first part of eq. ([78|) as 

nV^ 1 1 



(n/2) 

7+ 



^(_"/2) 



[1 



t = 



no 



(3/2) 

7+ 



V = 



1_ 

7+ ■ 
< V <l. 




(127) 
(128) 

(129) 
(130) 



Note at this point, that the behavior of p„ for large values of n is controlled by rtC„, a quantity 
that can have a kind of critical behavior in the present model. Namely, the large n behavior of 
nCn depends on the ratio of i = no/7^^^\ since for large values of n, we always have w" = 
(7-/7+)^"^^^ « 1- One may introduce a critical value of tiq, indicated as 

(3/2) 

(131) 



(3/2) 

7+ 



1 + x + Vl + 2a; 



and one may observe at this point that if no < tt-c one has limn^oonCn = 0, if no > nc, one obtains 
limn^oonCn = 00 and finally lim,n—,oanCn = 1 if no = ric- We shall return to the interpretation of 
this critical value of no later on. We shall see that the quantities C„ correspond to the combinants 
of the multiplicity distribution and no = ric critical value corresponds to the divergence of the mean 
boson multiplicity. 

Let us now proceed further with the solution of the model equations. The only remaining unsolved 
recurrence is eq. (^9|) which can be formally solved in a general manner, as follows. Let us introduce 
the generating function of the lj„-s and as Gui{z) and let us define an auxiliary F{z) as follows: 



F{z) 



U!„Z 



00 

^(n+l)C„+iz", 



and observe that 



F{z)G^{z) 



dtF{t)G^{t) 



E 

00 ^ 

E- 



Y.{k + l)Gk 



.fe=0 



n=l 



and this last equation can be re-written with the help of eq. ( |79| ) as 

dtF{t)G^{t) = G^(z)-1, 



F{z)G^{z) = 
GM) = 



az 

G^(0) exp 



dtF{t)) , 



(132) 
(133) 

(134) 
(135) 

(136) 
(137) 
(138) 



since the expansion coefficients of F{z)-s are known. The last step in the formal solution of eq. ( |79| ) 
is the introduction of the generating function for the probability distribution p„ as 



G(z) 



n=0 



(139) 
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which can be expressed with the generating function of the w„ distribution with the help of eq. ( fz^ ) 
as 

GM ^ (140) 

which yields the formal solution for the generating function of the probability distribution of p„ as 
follows: 

Giz) = exp|^f;C„(z"-l)^ , (141) 
where the quantities C„-s are the so called combinants of the probability distribution of p„ and in 



our case their explicit form is known for any set of model parameters, as given by eqs. (127,115,102 
|88| , p9| ). Note that this generating function depends on two variables, the parameter uq that controls 
the phase space density and the parameter x = R^jjcr^, eq. (102), that yields the available dimen 



sionless phase-space volume. Alternatively, one may introduce another set of variables by defining 
t = no/ric and v — vT-Tt+j < w < 1, to obtain a more transparent form of the probability 
generating function: 

Giz;t,v) = exp(^^^-^^— ^(z"-l)j. (142) 

We find that this multiplicity generating function is not corresponding to known discrete probability 
generating functions after inspecting a few standard textbooks on statistics that include a large 
number of probability generating functions, e.g. ref p^ . Thus we may assume that we have found a 
new type of probability generating function. We shall investigate its properties not only numerically 
but also analytically in certain simple limiting cases. Especially, we show that there is a critical 
value for hq at which a wave-packet version of Bose-Einstein condensation occurs that influences 
drastically the pn distribution too. We shall also show that this distribution reduces to known type 
of distributions in certain limiting cases. 

The combinants were introduced to statistics in refs. pTf and some of their properties were discussed 
recently in refs. p8[-pl|. The auxiliary quantity a;„ is the probability ratio 

l^n = Pn/PO- (143) 

(Note that in our model po = exp(— C„) 7^ ). The combinants up to a given order n can be 
expressed as a combination of the first n probability ratios as 

1 

Ca = UJo ^nCnUJq^n, (144) 



y„ — li/ 



n=l 



hence their name. It is outside the scope of the present paper to study in detail the properties 
of the combinants and their relationships to other characteristics of the multiplicity distributions, 
like cumulants and factorial moments, and scaling laws of count probabilities. For a more detailed 
discussion on these general topics, we recommend refs. p8|-[^. 

Let us note, however, a general property of a multiplicity distribution that is given in terms of its 
combinants. Any probability generating function that is given in terms of its combinants that are 
non-negative can be re- written as a convoluted Poisson distribution, or compound Poisson distribu- 
tion if Po > as 

G{z) ^ exp(C{H{z) - 1)), (145) 

where C = X^^^i can be interpreted as the mean multiplicity of the Poisson-distributed clusters 
(or clans). For completeness, we note that the probability generating function for a simple Poisson 
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distribution is Gp{z) — exp[n(z — 1)]. The particles within a single cluster are distributed according 
to the probability distribution 

^'a^ = (146) 

ji=i 

and there is always at least one particle in any cluster i.e. Pn=o ~ 0. 

Another interesting property of this multiplicity distribution is that the generating function can be 
written as 



Giz)^l[expiCniz''~l)) (147) 

n=l 

i.e. as an infinite product of generating functions of Poisson distributions of particle singlets, dou- 
blets, triplets etc, having a mean of Ci, C2, C3 etc respectively. The corresponding multiplicity 
distribution can therefore be expressed in terms of Poisson distributions of particle rt-tuples with 
means of Cn in a multiple convoluted manner |2^, ^ . This is a general property of any distribution 
that is given in terms of its combinants. 

Due to this property, the mean multiplicity can be expressed in terms of combinants as 



V"/ = ^ - 

n i—1 

We see that the finiteness of the mean multiplicity is related to the vanishing values of iCi for large 
i. This limit, on the other hand, can only be reached if uq < Uc- Thus one finds that (n) < cxd if 
no < Uc and (n) = 00 if np > ric. Thus Uc can be interpreted as a critical value for the parameter 
TiQ. We shall see that divergence of the mean multiplicity (n) is related to condensation of the 
wave-packet modes with the highest multiplicities. 

In general, the factorial cumulant moments, fq-s of the probability distribution can be expressed 
with the help of the combinants in a relatively simple and straightforward manner similarly to the 
results of refs. js^, |3^ : 

= ^ lnG(z)U.i = ^ * • - 1) • ... • (z - g + 1) • a (149) 

that are nothing else than the factorial moments within a single cluster multiplied with the average 
cluster multiplicity ||3C[| . 

Note, however, that the multiplicity distribution we find can not only be considered as an infi- 
nite convolution of Poisson distribution of particle singlets, pairs triplets etc but can be re-written 
also as an infinite convolution of Bose-Einstein (or Negative Binomial) distributions with coupled 
parameters. 

In the exponent of G{z) one may apply a negative binomial expansion of the terms 

1 ^ 1 ik + l){k + 2) f 

and with the help of this transformation, the generating function can be equivalently written as 
another infinite product: 



-I— 1- (fc + l)(fc + 2) 

G{z) = l[G,{z)^^, (151) 

k=0 
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This expression indicates that the probabihty distribution p„ can also be considered as an infinite 
convolution of generalized Bose-Einstein distributions. The mean multiplicities of the single modes 
increase in a geometrical series. For any single mode of this convolution, the probability distribution 
can be written as 



{nk + 1 )("+!) ' 

and the mean multiplicity for that given mode can alternatively be written as 



no 



no 



(154) 



(155) 



Thus the probability distribution can be considered both as a superposition of Poisson distributed 
independent clusters as well as an infinite convolution of Bose-Einstein distributions with coupled 
mean multiplicities. 

It is especially interesting to note, that in the dense Bose-gas limit, the fc = term dominates 
eq. (151), while in the rare Bose-gas limit each factor contributes with a similar weight from the 
infinite product in this equation. In contrast, one may consider the Poisson cluster decomposition 
as given by eqs. (145, 141) and in this type of decomposition, the first factor contributes to the rare 
gas limiting case but in the dense gas limiting case all factors are important. 

Let us make a few further remarks about the structure of the probability distribution p„ in this 
class of models. Although this multiplicity distribution is already given in terms of its generating 
function by cq. (142), and in terms of its combinants, Cn as given by eq. (127), a more explicit form 
for the multiplicities Pn can also be given. In general, if the combinants are known, the probability 
ratios a;„ = Pn/po can always be expressed as 



k=l 



ni— 1 n2 — 1 

E E 

n2—k — l n-^—k—2 



E^ 

nfc = l 



n— n2 C'n2— "3 



(156) 



Although this equation is seemingly more complicated than the recurrence given for ujn in eq. ([79|), 
it is still useful since the explicit form for C„ can be inserted from eq. (127) to find that 



ni— 1 712 — 1 



^"^E^^ E E -E 



fc=i 



712— k—l ns — fc— 2 



fc-l 

i=i 



-1/3 



where one can see that the t = no/uc dependence of w„ cx p„ is simple. This can be utilized to 
simplify the general evaluation of the LUn probability ratios as 



Fo{v) = 1 



Fn{v) 



1 " 

-E' 



n ^ — ' /i 



(158) 
(159) 

(160) 



This is the simplest form that we can find for probability distribution in the general case (the 
probabilities are given by eq. ( [7^ ) ). 

We investigate the probability generating function G{z) as well as the multiplicities and the particle 
spectra in more details in the rare and in the dense gas limiting cases in the next sub-sections. 



5.1 Analytic Results for Rare Bose Gas 

In this section we consider large source sizes or large effective temperatures, i.e. we examine the 
combinants and the generating function of the probability distribution in the a; >> 1 limiting case. 
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Since the multiplicity distribution was obtained analytically in terms of the combinants, let us 
evaluate these to leading order in 1/x « 1. One obtains that 



7+ 

7- 

7? -7- 



X 

n \ — 



— -n \ — 



2n 



(161) 

(162) 
(163) 

(164) 



These relationships can be substituted to the solution of the recurrence relations, eqs. (125-126) to 
obtain that 



1 



1 



gn 





)(3/2) ^3/2 


X 


^eff 




n 


no ( 


2\ 1(1-") 
xJ 



1(1-") 



(165) 
(166) 
(167) 



Note that the determination of a„ is a bit tricky, since a finite sub-leading part remains that is difficult 
to determine as the leading orders cancel. However, there is an important constraint between C„ 
and a„ as given by eq. (|7q). This constraint is satisfied by 



2(7^ 



2n 



Inserting these to the definition ( [77| ) one obtains that 
G„(ki,k2) = 







(riTTCTji) 



2 \3/2 



exp 



2ar 



-(k?+k^) 



2n 



(ki 



(168) 



(169) 



Both the inclusive and the exclusive momentum distributions will be built up from the auxiliary 
quantities Gn(ki, k2), and the leading order behavior is given by Gi(ki, k.2) in this rare gas limiting 
case, with first order corrections from G2(ki, k2). 

Eq. (167) indicates, that the leading order result for the combinants in the a; >> 1 limiting case is 



G„ 



and the first sub-leading correction is given by 

Cn = <5i,„no -I- <52,i 



nl 



2{2x) 



3/2 



(170) 



(171) 



Thus, to leading order, the probability distribution can be considered as a Poisson distribution 
of singlets with a sub-leading correction that yields a convolution of Poisson-distributed doublets. 
Similarly to this, the Poisson-doublets will modify the single-particle momentum distribution as 
well as the two-particle correlation function. We consider these modifications in two manner: i) hy 
evaluating the leading order corrections only and ii) by summing up certain sub-leading corrections. 
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5.1.1 Very rare gas 

In the rare gas limiting case, the leading order Poisson distribution can be utilized to solve the model 
completely as follows: 

n! 

n" 

Pn = — Texp(-no), (173) 



^^(^^ = ^'^^"^^ = (2.^Tt,)(3/^) (-^) ' (^^^^ 
P2(ki,k2) = p(^)(ki)Pf)(k2) [l + exp(-ii2^/(ki-k2)2)], (177) 
C2(ki,k2) = l + exp(-ii2^/(ki-k2)2). (178) 

(179) 

Thus the very rare gas limiting case is very simple and solvable completely in an analytical manner. 
The multiplicity distribution is a Poisson distribution, the mean multiplicity coinciding with no- 
The stimulated emission does not influence the probability distribution for the "unsymmetrized" 
Po,n that we assumed to have a Poisson form. There is only a normalization factor that appears as a 
modification to the single-particle momentum distributions, and, apart from an overall normalization 
factor Afc the two-particle inclusive correlation function takes a Gaussian form with an eflEective radius 
parameter Reff- For large mean multiplicities, the overall normalization factor approaches unity 
since lim„o^oo AAc = 1. 

To our best knowledge, no similar evaluation of the two-particle inclusive correlations from the 

recurrence relations of S. Pratt was performed before, only exclusive correlation functions were 
evaluated numerically Thus, this is the first time as far as we know that the exclusive momentum 
distributions are summed up and the inclusive particle spectra and the inclusive correlations are 
determined from the plane-wave model, too. Our method seems to be powerful enough to do certain 
analytical calculations even for systems where the numerical evaluation of the inclusive distributions 
were too tedious to perform earlier. Let us also stress that these results correspond to a multi-particle 
wave-packet system. 

One may argue that the solution presented above is rather trivial and that not much is learned about 
the nature of the symmetrization effects here since they cancel from the leading order results in the 
rare Bose-gas limiting case. But the above simple results became very useful when one determines 
the next to leading order results in the rare gas limit, and it is necessary to know what is the reference 
point and what kind of changes happen if the leading order contribution from the symmetrization 
is taken into account. 

5.1.2 The rare gas limiting case 

Let us evaluate now the leading order corrections to the Poisson limiting case. This corresponds to 
keeping the first order corrections in 1/x. The probability generating function reads as 



G{z) = exp {n„{z - 1) + C2{z^ - 1)) ^ exp (no(z - 1)) (^1 + -^-^L(^2 _ 
which yields the following multiplicity distribution: 

"■0 / \ 

Pn = — r exp(-no) 



n(n — 1) — Hq 
2(2x)(3/2) 



(180) 



(181) 
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We see that the leading order Poisson multipHcity distribution is modified in such a manner that 
the corrected multiphcity distribution is depleted in the small n region, where n{n — 1) < ?t.q and 
enhanced in the high multiplicity region, where n{n — 

Pn , n{n - 1) 



^0' 



= 1 



(182) 



Po,n 2(2x)(3/2) 

This shift to large multiplicities is more enhanced for smaller values of x i.e. for higher phase-space 
density. The mean multiplicity can be calculated as 

^ ^2 



dz 



no 



(183) 



(2x)(3/2) 

which increases in this leading order calculation quadratically with increasing values of no due to 
the influence of Poisson-doublets. Note also the strong rise of the mean multiplicity with decreasing 
values of x i.e. with increasing phase-space density of the wave-packets. 

Evaluation of similar leading order corrections to the spectra and the correlations is straight-forward. 
The n-particle exclusive and inclusive number-distributions have simpler form than the correspond- 
ing probabilities, hence we indicate here these quantities: 

k2 



7v{")(k) 



A^i(k) = 



(7r4)(3/2) 

n{n — 1) I 
(2a;)(3/2) ] 

(7r4)(3/2) 



exp 



1 (3/2) 



ttctt 



exp 
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(3/2) 
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exp ( 













exp 



2 



exp 



(184) 



(185) 



(186) 



where the mean multi- 



(^4a;)(3/2) 

Note that these quantities are properly normalized according to eqs. (^ 

plicity is given by eq. ( |183| ). 

The interpretation of these results is straight-forward: i) The single-particle inclusive momentum 
distribution is enhanced at low momentum. 



The enhancement is proportional to 



;:372 



which is the mean number density 



of pairs in the phase-space volume that is available for a single particle in the rare gas limit. 
ii) The single-particle exclusive momentum distribution of rank n is enhanced at low momentum, 
the enhancement being proportional to the number of pairs, n{n~ 1). This predicts a characteristic 
dependence in the exclusive momentum distribution as a function of n, that is due to the Bose- 
Einstein symmetrization. Such dependences could be studied in the event by event analysis to be 
performed by the NA49 Collaboration, for example. 

The evaluation of the Bose-Einstein correlation functions in this limiting case is rather involved and 
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Let us turn our attention to the opposite limiting case, the ultra-dense Bose-gas of wave-packets. 



5.2 Analytic Results for Highly Condensed Bose Gas 

This limiting case corresponds either to a small effective radius parameter or a low effective tempera- 
ture of the multi-boson wave-packet system. Formally, it corresponds to the x = 2mTef fR^^ j << 1 
limiting case. In this case, one can determine the leading order expression for the combinants of the 
probability distribution as 

7+ = l + x-^0{x^), (187) 
7_ = + 0(2:2), (^gg) 



nc{x) ~ (1 + a;) 



3/2 



no 

Ur 



1 - 

3x 



nl + 0{x''). 



(189) 
(190) 
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Thus the probability generating function can be written in the condensed Bose-gas limiting case as 



- (1-^(1-.)) ^'''^ 
n = (192) 

which corresponds to a Bose- Einstein distribution in the a; << 1 limiting case with a mean multi- 
plicity that diverges if no reaches the value of — (1 + 3a;/2). In this a; << 1 limiting case, the 
multiplicity distribution shall be given by a Poisson distribution, and the model can be completely 
solved analytically as follows: 

n" 

Pn = ^^^TTJTT' (193) 



(n+ 1) 



nj {n + iy 



(194) 



GuM = (_^^^„+,G,(p,q) = (^ + i).+i^;;^|:)(^-p(^-^^j' (195) 

Pi(k) = p(")(k) = — — L_expf-— (197) 



(27rmTe//)(3/2) V 2mTe// 

P2(kl,k2) = Pi(ki)Pi(k2), (198) 

C72(ki,k2) = 1. (199) 

Note that this solution is obtained by summing up leading order contributions in the small x region. 
In this limiting case, corresponding to a; = 2mReffTeff << 1, essentially all the wave- packets remain 
condensed in the same wave-packet state as reflected by the shape of the single-particle inclusive 
spectrum and by the vanishing enhancement of the inclusive Bose-Einstein correlations. The mean 

3x 

multiplicity diverges as no approaches nc — 1 -I — "^^^^ corresponds to a very small wave-packet 

system where stimulated emission can dominate the particle emission even for very small values of 
the parameter no. The multiplicity distribution is a Bose-Einstein distribution, the parameter no 
can be interpreted as a measure of Pn/Pn+i < 1 in units of Uc, the inclusive and exclusive single- 
particle distributions coincide and there is no enhancement in the two-particle inclusive correlation 
function. Due to this reason, the solution corresponds to a completely coherent behavior, when 
all the particles are emitted in the same wave-packet state. It is interesting to observe, that the 
coherent behavior in this picture (defined by the vanishing enhancement in the two-particle inclusive 
correlation function) appears simultaneously with a thermalized Bose-Einstein distribution of the 
multiplicity distribution for ki ^ k2. Thus the effective coherence appears due to the negligibly 
small effective source-size Reff — > 0, or a very small efi'ective temperature T^ff 0. In these 
limiting cases, the model corresponds to a thermalized source prepared in such a manner that all 
the wave-packets are created in the same state and their overlap is maximal. This is the reason 
why thermal behavior and coherence can appear simultaneously in this picture. It is interesting 
to note that the single-particle inclusive and exclusive momentum distributions coincide and are 

2 

described by a Boltzmann distribution with an effective temperature T^ff =T+ that picks up a 

contribution from the temperature of the source and from the momentum- width of the wave-packets. 
This shape of the momentum distribution is forced by the Gaussian ansatz for the momentum-space 
representation of the wave-packets and the Gaussian distribution of the mean momentum of the 
packets. 

In this and the previous sub-sections, analytical solutions to the multi-particle wave-packet system 
are presented. These analytic solutions yield new result not only in the study of multi-particle wave- 
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packet problems, but also of the plane-wave problem that corresponds to the wave-packet problem, 
obtainable by the substitution Reff R and Tgff — s- T. 



5.3 Other Results for the Critical Multiplicity 

The critical behavior, related to the divergence of the mean multiplicity in the original plane-wave 
system was discussed first in ref. Q]. In this pioneering work, S. Pratt published a formula for the 
onset of the condensation, determining the critical number of the unsymmetrized pions from the 
inspection of certain ring-diagram. Although Figure 1 of ref. § was evaluated with the correct 
expression p3. 



(3/2) 



(200) 



this expression appeared in the paper in a misprinted manner, eq. (9) of ref. Q, where r]c = ric and 
A = mT. Note that the steps to derive eq. (9) were not given in ref. [Q and the misprinted formula 
depends also on po = /{2m). Our result on the critical density for the condensation we define can 
be written as 

(3/2) 



71c = 



1 + 2mTeffR, 



l + 4mTeffR 



(201) 



which corresponds to the correct formula of eq. (200) when going to the wave-packet system from the 
plane- wave system with the help of the replacements rjc ^ A — mT —> mT^ff and R — *■ Reff- 
On Figure (|l|) we indicate the critical pion multiplicity as a function of wave-packet size. 
It is important to discuss how the critical value of the parameter tlq can be determined. We were 
interested in finding that limiting case, when the explosion of the (symmetrized) Pn probabilities 
just happens. This we call the Bose-Einstcin condensation point and we define it physically as the 
curve in the n, i?e//, T'e// parameter space along which the decrease of the (unsymmetrized ) pn^ 
probabilities is just compensated by the stimulated emission of bosons. With other words, this is the 
point where a strong stimulated emission is build-up, which corresponds to the onset of the lasing 
mode. 

Thus the condensation point is defined as 



npn{nc 



{n + l)pn+i{nc) 



(202) 



which is to be satisfied in the large n limit only. It is clear from this definition, that when this 
condition is satisfied, the mean multiplicity has to diverge. 

Let us first determine the critical boson multiplicity for the onset of the Bose-Einstein condensation 
of the wave-packets into the same wave-packet state. From npn = {n + l)pn+i in the large n limit it 
follows that nUn = (n -\- l)ti;„+i. From this it follows that nC„ = {n + 1)C„ in the large n limiting 
case. Comparing this to eq. (115), one finds that at the onset of Bose-Einstein condensation, as 
noted before, happens at 



no 



1 -I- X -I- Vl + 2x 



(3/2) 



(203) 



where x — R^^j^a^ — 2mTef fR^j j is a dimensionless measure of the size of phase-space cell occupied 
on average by a single quanta. This expression is identical to eq. ( |20l| ). One may determine this 
critical Uc parameter of the distribution as a function of {Reff, Teff} or as a function of {R, T, cr}. 



Thus the parameter no can be compared to its critical value Ue and if no < nc, lim„ 



but 



above the critical multiplicity no > lim„^oop„ = oo i.e. stimulated emission over-compensates 
the decrease of the unsymmetrized emission probabilities and a coherent, laser-like behavior occurs. 
Thus the ratio no / n^ controls the competition between the stimulated emission and the decrease of 
the production probabilities of the unsymmetrized emission, and n^ can be interpreted as a critical 
pion density at which the condensation of wave-packets starts to appear. 
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6 Numerical Results for Mult i- Particle Wave-Packets 



In this section we present some numerical results about the effects of multi-particle symmetrizations 
on the multiplicity distribution, on the single-particle exclusive momentum distribution and on the 
two-particle exclusive correlation functions. Although these results are more readily interpretable 
when one can rely on the analytical insights gained in the previous section, a qualitative understand- 
ing of the multi-particle symmetrization effects on these observables can be obtained from a direct 
inspection of these figures. Let us investigate first the modification of the multiplicity distribution, 
Pn, due to the overlap of the wave-packets. 

In Figs. H and || the multiplicity distribution, eq.(60) is plotted for various wave packet sizes. The 
"multiplicity parameter", no, was fixed to ng = 600. These Figures clearly demonstrate, that the 
distribution function approaches a divergent region for wave packet sizes in the vicinity of ax = '^fm 
for this set of parameters. This is due to the fact that the overlap between the wave-packets is 
controlled by the wave-packet size and if the overlap reaches a critical magnitude, stimulated emission 
of wave-packets in similar wave-packet states leads to a multiplicity distribution that diverges for 
large values of n. As the size of the overlap of the multi-particle wave-packets is increased, the 
multiplicity distribution becomes more and more different from the Poissonian limit, and it is shifted 
to larger multiplicities in accordance with our leading order analytical result given by eq. ( |181| ). 
Investigating the single-particle exclusive momentum distributions, our analytical results about the 
development of a low-momentum enhancement is confirmed also numerically, as can be seen on 
Figures ^ and |^. It is also worth mentioning that the modification of the slope parameter of the 
high momentum part of the distribution is significant only when the spatial width of the wave- 
packets is small. This happens since the effective slope parameter, T^ff is given by eq. and the 
contribution of the wave-packets to this quantity is large only if the momentum- width of the packets 
is large i.e. their spatial width is small. 

The effect of the wave packet size on the two-particle exclusive correlation functions is shown on 
Figs. H and 0. 

The virtual size of the source, i?^, is assumed to be Ry = l/|Ak|]^/2> where |Ak|x/2 is the momentum 
difference, at which the correlation C(Ak) — 1 drops to half of its intercept value, C(0) — 1 .Our first 
expectation would be, that Ry is a monotone increasing function of (Jx- However, as Fig. |^ shows, 
this expectation is not fulfilled. This effect is even more pronounced in Fig. |^. The explanation of 
this behavior can be given as follows. From Figs. 3 and 4. we have seen, that the system approaches 
the formation of a condensate at ax ~ 4.0 fm. At this size, the overlap between the wave packets 
within the "virtual volume" is maximal. This case corresponds to the smallest "virtual" radius 
parameter of the exclusive correlation function and to the smallest value of its intercept parameter 
A = C(k, k) - 1. 

It is also clear, that the condensation effects the low momentum pions. Thus as we changed the 
mean pion momentum from 100 MeV to 50 MeV, we included more and more from the condensed 
part, which is responsible for the strange behavior. 

Note that the "virtual" radius parameter first decreases than increases with increasing values of the 
wave-packet size ax ■ Similarly, the intercept parameter X = C2 (k, k) first decreases then it increases 
back to its conventional value of 1, for increasing size of the wave-packets. As we have seen in the 
ultra-dense limiting case, A = is also possible for high densities. 

7 Summary, Conclusions and Outlook 

In this paper a consequent quantum mechanical description of multi-boson systems is presented, 
using properly normalized projector operators to overlapping multi-particle wave-packet states for 
bosons. Hundreds of overlapping wave packets lead to a difficult problem. We found, however, a 
possible configuration, for which an exact solution is obtained. In our description no phase averaging 
is used. The mathematical problem imposed by the large number of source points is dealt with an 
algebraic procedure, similar to that invented by S.Pratt in ref. |l[]. The effects arising from the 
multi-particle symmetrization and from the finite width of the wave packet, ax, are shown for a 
system that has a radius R and a freeze-out temperature T similar to a fireball that could be formed 
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at mid-rapidity in Pb + Pb collisions at CERN SPS. The effects that depend on the wave-packet 
size are as follows: 

1) . The critical density of pions, at which the condensation appears, is a function of , and one 
finds that a pion-laser is formed at around = 4: fm, when assuming tiq — 600, i?, = 11 fm and 
T = 120 MeV. 

2) . The multiplicity distribution shifts toward higher and higher multiplicities as is increased 
towards 4 fm. In a range of wave-packet widths that include 4 fm, the multiplicity distribution ex- 
plodes, i.e. we approach the condensation point. As is further increased, the average multiplicity 
decreases back again to normal values. 

3) . The single particle momentum distribution shows a very large inverse slope parameter for small 
cr^y say (Tj. = 1 fm. This is a trivial effect: in the Fourier transform of a narrow wave packet there 
is a large contribution from high momentum components, which shows up in the single particle 
momentum distribution. Further, as approaches the 4 fm value, the low momentum part in the 
spectrum becomes more and more enhanced. This is understandable, if we keep in mind, that the 
Bose condensation is approached first for the low momentum part of the spectrum. As is further 
increased, the spectrum approaches again a Boltzmann distribution. This can be understood, if 
we take into account, that for larger the virtual volume increases and thus the density of pions 
within this volume decreases. 

4 ) The two particle correlation function shows an interesting behavior. Two features are to be 
emphasized: i) The virtual radius first decreases and later increases with increasing ax- This effect 
is larger as the mean momentum of the two pions decreases, ii ) The intercept value, C( | ki — k2 1 — 0) , 
decreases as is increased from 1 fm to 4 fm, and increases again as is further increased. These 
effects are caused by the increase and decrease of pion density within the " virtual volume" , which is 
influenced by the distance to onset of the divergence of the mean multiplicity, i.e. the onset of the 
laser (or condensation) mode. 

Finally, we may conclude, that we have found interesting effects for the investigated system. The 
construction of a clear connection of these effects to the onset of Bose-Einstein condensation as it is 
known in statistical physics, is to be made. Further, the problem, whether our approach is a useful 
one for the Bose-Einstein condensation of atomic systems, remains also for future investigations. 
Such studies may be promising since the atoms have to be trapped with the help of a magnetic trap 
to cool them below the critical temperature. The usual theoretical description uses the non-linear 
Schrodinger equation where a potential is created from \4>\'^ ■ Because of the trap, plane wave states 
cannot be utilized. Thus the wave-packet description, presented in this paper may have a relevance 
also in the field of Bose-Einstein condensation of atomic vapors. 
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Appendix A 



Here we collect some formulas for another type of factorization of the n-particlc density matrix. 
For the n-particle density matrix we also may assume, that it is factorizablc in a simple way: 

n 

p„(ai, ...,q;„) = Jl p(ai) (204) 

4=1 

This choice of the density matrix is one of the possibilities. In the body of the paper we study an 
other case where the integrations can be carried out much easier, referred to as "model A" in the 
forthcoming. In this appendix, we explore the consequences of "model B" , where the integrals are 
much more difficult and the calculation can be performed only in some approximate manner. 
Although in this general case one obtains quite complicated integrals when evaluating the momentum 
distribution, because of the overlapping of the wave packets, we can still proceed a little further 
with these expressions. Both in case of model A and B, these integrals for the m-pion momentum 
distribution for fixed multiplicity n contain the factor given by Eq. ( ^ ) . 
In case of model B, the two-particle momentum distribution is given as 



/n 
darn P{ani) 



w*(ki,Q;j)u)*(k2,Q;j) w(ki,Q;crJ w(k2,ao.J 7i,cr, 



(205) 



n 7fe,-. 

cr(") fe=l 



while in case model A, the integrals become somewhat simpler due to the cancelation of the extra 
factor in the density matrix with the denominator of the expectation value in eq. ^ 
For model B, the single particle distribution reads as 

n n 

^ ^w*(ki,a4)-!«(ki,aaj H ^'•'^^ 
<i(ki) = [ f[da^p{a,n)^^^ , (206) 

En 7^-.. 

o-(") fc=l 

In case of model B, the integrals defining the spectrum and the two-particle distributions can also be 
reduced to a superposition of integrals, which contain only Gaussian factors, using an expansion of 
the denominator in the integrand. When all -s are equal, the denominator reaches its maximum 
with maximum value of n\. Thus one can expand it in the absolutely convergent series: 



1 1 



EtxC) riLi 7fc,^. nl{l-x) 
where we introduced the quantity < a; < 1 with the definition 



-. oo 

-Y^.' (207) 



(208) 



This expansion is absolutely convergent in terms of x in its domain. Since contains linear 
combinations of different powers of 7 factors, which are themselves Gaussian factors, their powers 
are also Gaussian and since all the remaining factors in the integrals were Gaussian, we are left with 
Gaussian integrals. However, these results are rather complicated even in the n — 2 case and we do 
not include them in the present work. Due to the greater analytical simplicity of model A, we have 
explicitly evaluated integrals only for this latter case, as was given in the body of the paper. 
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Figure Captions 



Fig. 1. Solid line stands for the critical pion multiplicity distribution, ric as a function of the wave 
packet sizes a^, as indicated by solid line. The other parameters have the following values: T — 120 
MeV, i? = 11 fm, no = 600. 

Fig. 2. Pion multiplicity distribution for wave packet sizes a^; = 1.0, 2.0 and 2.5 fm, as indicated by 
solid, dashed and dotted line, respectively. The other parameters have the following values: T = 120 
MeV, R=n fm, no = 600. 

Fig. 3. Pion multiplicity distribution for wave packet sizes ax = 4.0, 7.0 and 10.0 fm, as indicated 
by solid, dashed and dotted lines, respectively. The other parameters have the following values: 
T — 120 MeV, R — 11 fm, uq — 600. At a critical size of the wave-packets, the overlap between 
the packets will be sufficiently large to start an "alavanche" of induced emissions, characterized by 
a Pn distribution which increases with increasing values of n. (Such distributions can be normalized 
if they are truncated at large values of n). As the spatial width of the wave-packets is increased 
further, they become more and more similar to plane-waves in momentum space and the overlap of 
the wave-functions is decreased. The multiplicity distribution is shifted back towards the Poissonian 
limit. 

Fig. 4. Exclusive single particle momentum distribution, p|^°"-'(k) for a sub-set of events with 
a fixed multiplicity of n = 600. Eq. jl^ ) is plotted for wave packet sizes = 1.0,2.0 and 2.5 
fm, as indicated by solid, dashed and dotted lines, respectively. The other parameters have the 
following values: T — 120 MeV, i? = 11 fm, no — 600. The x-axis is scaled with k^. Note that the 
effective slope parameter changes with the variation of the wave-packet size. As the critical overlap 
is approached, the low momentum peak becomes more and more pronounced in the spectrum. 

Fig. 5. Exclusive single particle momentum distribution, P^^^^\'k) as given by eq.dl^) for wave 
packet sizes a = 4.0, 7.0 and 10.0 fm, as indicated by solid, dashed and dotted lines, respectively. 
The other parameters have the following values: T = 120 MeV, i? = 11 fm,no = 600. The x-axis 
is scaled with k^. The slope parameter at high momentum is hardly changed when ax is varied in 
this range, but the low momentum enhancement is strong only if the overlap of the wave-packets 
approaches the critical value. 

Fig. 6. Two-particle exclusive correlation function, as given by eq. (p7|) , for wave packet sizes 
ax = 1.0, 2.5 and 10.0 fm is plotted with solid, dashed and dotted lines, respectively. The mean 
momentum is fixed to |ki -I- k2|/2 — 100 MeV. The relative momentum is parallel to the mean 
momentum, corresponding to the out direction for spherically symmetric systems. The x-axes is 
scaled with the relative momentum, |ki — k2| = Ako„t. The number of pions was fixed to n = 700. 
The "virtual radii" are Rv{ax = 1.0 fm) = 13.2 fm, R.^{ax = 2.0 fm ) = 11.6 fm, Ry{ax = 2.5 fm 
) 15.2 fm. 



Fig. 7. Two-particle exclusive correlation function, C2^°°^(ki, k2) is plotted, as given by eq.(17), 
for wave packet sizes ax = 1-0, 2.5, 4.0 and 15.0 fm with solid, dashed, dotted and dense-dotted 
lines, respectively. The mean momentum is fixed to |ki -I- k2|/2 = 50 MeV. The two momenta are 
parallel (out component). The x-axis is scaled with the momentum difference, |ki — k2|. The actual 
number of pions was set to n = 800. The virtual radii are Rv{ax — 1-0 fm ) = 12.33 fm, Rv{ax = 2.5 
fm ) = 8.57 fm, i?^(a^ = 4.0 fm ) = 8.22 fm, Ry{ax = 15.0 fm ) = 17.4 fm. 



31 




Figure 1: 




Figure 2: 



32 



0.04 



0.03 



0.02 



0.01 



= 10 fm 



300 400 500 



600 
n 



a, = 4 fm 




700 800 900 



Figure 3: 



le-07 



le-08 



le-09 



a, = 2.5 fm 



= 1 fm 



= 2 fm 



10000 20000 30000 40000 50000 60000 70000 80000 

/MeV^/ 



Figure 4: 



33 



le-07 



le-08 



le-09 





= 4 fm 




~ a^ = 7fm 
















= 10 fm 



10000 20000 
r [MeVT 



30000 



40000 



Figure 5: 



2.2 




0.8 



10 20 30 40 

ki - k2 [MeV] 



50 



60 



Figure 6: 



34 




kj - [MeV] 
Figure 7: 



35 



